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Abstract 

We compute the complete third-order contributions to the coefficient functions for the longitudinal 
structure function Ft, thus completing the next-to-next-to-leading order (NNLO) description of 
unpolarized electromagnetic deep-inelastic scattering in massless perturbative QCD. Our exact 
results agree with determinations of low even-integer Mellin moments and of the leading small-x 
terms in the flavour-singlet sector. In this letter we present compact and accurate parametrizations 
of the results and illustrate the numerical impact of the NNLO corrections. 



We have recently published the complete three-loop splitting functions p(^) for the evolution of 
unpolarized parton distributions of hadrons |ni,l2l| . Together with the second-order coefficient func- 
tions ll3l|4l|5l|6n71, these quantities form the next-to-next-to-leading order (NNLO) approximation 
of massless perturbative QCD for the structure functions Fi , F2 and Ft, in deep-inelastic scattering 
(DIS). The situation is different, however, for the longitudinal structure function Fi = F2 — 2xF\. 
Here the leading contribution to the coefficient functions is of first order in the strong coupling 
constant tts (instead of as in ^1,2,3)^ and consequently the scheme dependence of the splitting 
functions P^'^^ is compensated only by the third-order coefficient functions. The latter quantities 
are thus required for completing the NNLO description of the structure functions in DIS. 

At this point there is only very limited (and no entirely model independent) experimental infor- 
mation on Ft in the most interesting region of very small values of the Bjorken variable x ISllQl fTOl 
(see ref. [11 J for an overview over previous measurements at large x). Dedicated runs of HERA 
with lower proton energies would significantly improve this situation ifT^ [T3 1 and provide valu- 
able further constraints on the proton's gluon distribution at very small momentum fractions. The 
NNLO corrections are required to assess in which kinematical region of x and = —q^, with q 
being the momentum of the exchanged gauge boson, such constraints can be reliably obtained. 

As discussed in refs. ||I1I3, see also refs. lfT4lfT5lfT6ll . our determination of the NNLO splitting 
functions P^'^'^ has been performed via a three-loop Mellin-A^ space calculation of DIS in dimen- 
sional regularization with D = 4 — 2£. The splitting functions are then given by the coefficients of 
the 1/e poles, while the £^ results include the third-order coefficient functions. Hence by keeping 
the terms of order £^ throughout the computations, and by adding the second Lorentz projection of 
the hadronic tensor required to disentangle F2 and Ft, we are able to obtain the third-order coef- 
ficient functions as well. Our approach closely follows the calculation of the lowest even-integer 
moments in refs. ifTTlfTSlfTOl . Incidentally, the methods for obtaining all-A^ results at higher orders 
have been pioneered in the calculation of ref. |20| for the structure function Ft at second order. 

In this letter we present our jc- space results for the coefficient functions for Ft in electromag- 
netic DIS at three loops in a parametrized, but sufficiently accurate form. The lengthy exact for- 
mulae will be presented, together with the results for F2, in a forthcoming publication |21 1. 

Disregarding power corrections, the longitudinal structure function can be written as 

X^^Ft = Ct,ns®qm+ (e^) {Ct,q^qs+Ct,g^g) ■ (1) 

Here qi and g represent the number distributions of quarks and gluons, respectively, in the frac- 
tional hadron momentum, q^ stands for the flavour-singlet quark distribution, q^ = E;=i(^i + 
where denotes the number of effectively massless flavours, qns is the corresponding non-singlet 
combination. The average squared charge (= 5/18 for even Hj-) is represented by (e^), and de- 
notes the Mellin convolution. We write the perturbative expansion of the coefficient functions as 

C.,.(a„x)= i:^(^)"c2W. (2) 

Note that, as usual for Ft, the superscript of the coefficients on the right-hand- side represents 
the order in tts and not, as for the splitting functions, the 'm' in N™LO. Throughout this article 
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we identify the renormalization and factorization scales with the physical scale Q^. The results 
beyond the leading order (LO), n = 1 in eq. Q, are presented in the standard MS-scheme. 

For the convenience of the reader we first recall the known first and second order results. The 
scheme-independent LO coefficient functions for Fi read 

4!LW=4Qx, c£,W=0, c^ll{x)=%nfx{\-x), (3) 

with = 4/3 in QCD. Here and below, the singlet-quark coefficient function is decomposed into 

the non-singlet and a 'pure singlet' contribution, c^"q = c^"^^ + To a sufficient accuracy, the 
NLO {n = 2) longitudinal coefficient functions ll^ ffllTll can be written as 

^SsW - 128/9jcL^-46.50xLi-84.094LoLi- 37.338 + 89.53X 

33. 82^2 -f;cLo(32.90+ 18.41 Lo) - 128/9Lo- 0.012 5(jci) 
+ 16/27n^{6jcLi -12xLo-25x + 6} , (4) 

4^JsW - n/{(15.94-5.212x);c?Li + (0.421 + 1.520;c)Lg + 28.09xiLo 

- (2.370x"^- 19.27) } , (5) 

c[^^(j«;) ^ { (94.74 -49.20x)jciLf + 864.8;ciLi + 1161 xLiLo 

+ 60.06jcL^ + 39.66jciLo- 5.333 (jc"^-l)} , (6) 

where we have used the abbreviations 

xi = \ — X , Lq = \n X , Li — \n xi . (7) 

Eq. (Ill) improves upon the previous parametrization in ref. ll23ll . while eqs. Q and ^ have been 
directly taken over from ref. ll24il . The expressions © - as well as their convolutions with 
typical parton distributions, deviate from the exact results by one part in a thousand or less. 

Now we turn to our new third-order results. Again using the abbreviations (Q and inserting the 
numerical values for the colour factors, the non-singlet coefficient function can be parametrized as 

^SsW - 512/27 L|-177.40Lf + 650.6L?-2729Li- 2220.5 -7884x 

+ 4168;c2-(844.7Lo + 517.3Li)LoLi + (195.6Li-125.3)xiLf 
+ 208.3x4 - 1355.7 Lq - 7456/27 - 1280/81 + 0. 1 13 6(jci) 
+ ny{1024/81L^- 112.35 L? + 344.1Li +408.4-9.345^-919.3x2 
+ (239.7 + 20.63 Li) xiL? + (887.3 + 294.5 Lq - 59. 14Li)LoLi 
- 1792/8 1 xLI + 200.73 Lq + 64/3 + 0.006 5(jci ) } 

+ n} {3x11 + (6 - 25x)Li - 19 + (317/6 - 12^2) - 6xLoLi + 6xLi2(x) 
+ 9xLI - (6 - 50x)Lo} 64/81 

+ //rin^{(107.0 + 321.05x-54.62x2)xi-26.717 + 9.773Lo 

+ (363.8 + 68.32Lo)xLo- 320/81 L^(2 + Lo)}x , (8) 
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where the nj part is exact. The (non-//ii) fermionic contributions were already computed in 
ref. I.14J . The NNLO pure- singlet coefficient function for Fl can be approximated by 

4^JsW - {(1568/27 L^-3968/9L? + 5124Li) 4 + (2184Lo + 6059;ci)LoLi 

- (795.6 + 1036jc)4 - 143.6xiLo + 2848/9 - 1600/27 

- (885.53JC1 + 182.00 Lo)x"^xi} 

+ nj { (-32/9 Lf + 29.52 Li);cf + (35.18Lo + 73.06;ci)LoLi - 35.24xL^ 

- (14.16 - 69.84 - 69.41 jciLo - 128/9 + 40.239 jc^^jc?} 
+ flfl rif {(107.0 + 321 .05 X - 54.62 - 26.7 17 + 9.773 Lq 

+ (363. 8 + 68. 32Lo)jcLo-320/81L^(2 + Lo)}jc . (9) 

Finally the corresponding gluon coefficient function can be written as 

4^gW - '^/{(144L|-47024/27L^ + 6319L? + 53160Li)xi +72549 LoLi 
+ 88238 L^Li + (3709 - 335 14x - 9533;c2) xi + 66773 xLI 

- 1 1 17 Lo + 45.37 - 5360/27 - (2044.70jci + 409.506 Lq) x'^} 

+ n/{(32/3L^ - 1216/9L?-592.3Li + 151 1 ;cLi);ci + 31 1.3LoLi 
+ 14.24 L^Li + (577.3 -729.0x)jci +30.78 xL^ + 366.0Lo 
+ 1000/9L^ + 160/9L^ + 88.5037x"ijci} 

+ /Zfjn^ {(-0.0105 L^ + 1. 550 L? + 19.72jcLi-66.745x + 0.615jc2)xi 
+ 20/27 JcL^ + (280/8 1 + 2.260 jc)jc4 - ( 15 .40 - 2.201 x) jcL^ 
-(71.66-0.121x)jcLo} . (10) 

Eqs. © - (fTUI) involve the new charge factors (see refs. |Il7l[THll for the corresponding diagrams) 

/Zi7 = 3(e) , /Zfi = {e)^/{e^) , /Zff = flf, -fir, dD 

where (e^) stand for the average of the charge for the active quark flavours, (e^) = • 
The rational coefficients in eqs. dH) - (flOb are exact, as are (up to the truncation) the irrational 
coefficients of Lj in eq. (EJ) and of the 1 /x terms in eqs. Q - (flOt . The remaining coefficients have 
been fitted to the exact coefficient functions which we evaluated using a weight-five extension of 
the program f25] for the evaluation of the harmonic poly logarithms [26| (see ref. |27| for a new 
program without restrictions on the weight). Also the parametrizations dH) - dlOb deviate from the 
exact results by less than one part in a thousand, except for values close to zeros of Ci ^{x). 

Our results agree with the lowest six/seven even moments computed in refs. llTTlfTHlfT^ . The 
same holds for the sixteenth moments for both F2 and Fl computed in ref. 1.2^ as another check 
of our calculation. We also find agreement of the leading small-x terms x~^lnx of c^p^ and 
with the predictions obtained in ref. 1.291 from the small-.x; resummation. Moreover a threshold 
resummation has been derived |l30l|3T]| for the large-x logarithms in C^ns- We are however not 
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aware of any third-order predictions derived in this framework. Note that also at this order in 
tts the coefficient of the leading large-x logarithm of c^^ ( = 8 Cp in our notation - recall that 
we expand in terms of a^/ (47l) ) equals that of the leading -i-distribution of the corresponding 
coefficient function for F2 lfT8ll2l1 . We expect that this relation hold to all orders. Together with 
the prediction of the soft-gluon resummation for F2 (see, e.g., ref. [■321 ) this leads to the conjecture 

41W|_^ = ^["JWI^^ = ^^ln^"-^(l-.) + 0(ln2-3(l-.)) , (12) 

where, of course, the last term has to be replaced by O ( 1 — jc) for n=\. Eq. (fT2b will, e.g., be useful 
for approximate reconstructions of cj^q(jc) from a future four-loop generalization of ref. ifTTI . 

We now illustrate the size of the contributions Q - (flOb to eqs. (HI) and Q, for brevity focusing 
on the flavour singlet sector. The perturbative expansions of CL,q = CL,ns +Cl,ps and of C^g are 
shown in fig. 1 for four flavours and a typical value tts = 0.2 of the strong coupling. Especially for 
CL,g, both the second-order and the third-order contributions are rather large over almost the whole 
x-range. Most striking, however, is the behaviour of both CL,q and Ci.g at very small values of x. 
Here the anomalously small {xcf^ ~ x^) one-loop parts are negligible against the (negative) con- 
stant two-loop terms, which in turn are completely overwhelmed by the (positive) new three-loop 

(3) 

corrections xc\^ ^ ~ \nx-\- const . All we know beyond the third order are the leading small-jc terms 
derived in ref. ||29| . The resulting fourth-order (N^LO) corrections in turn exceed our three-loop 
contributions Q and (fTOb . Recall, however, that the leading low-x terms often substantially over- 
estimate the complete results at accessible values of x, as shown in fig. 1 for the third-order results, 
and that small-x information alone is usually insufficient for reliable estimates of the convolutions 
connecting the coefficient functions to observables as in eq. ([T]), see also refs. ID 13. Hence it is 
not possible to draw any conclusions about the behaviour of Fi beyond NNLO at this point. 

In figs. 2 and 3, CL.q and CL,g are convoluted with the sufficiently realistic model distributions 
xq,{x) = 0.6jc"°-3(1-jc)3-5(1+5.0x^-^) , 

xg{x) = 1.6jc"°-3(1 -jc)4-^(l -0.6x^-3) , (13) 

corresponding to a typical scale of about 30 GeV^. A comparison of figs. 1 and 2 clearly reveals the 
smoothening effect of the Mellin convolutions. In fact, under the chosen conditions, the (mostly 
positive) NNLO corrections to the flavour- singlet F^ amount to less than 20% for x < 0.3 and to less 
than 10% for 5 ■ 10^^ < x < 0.1. Note that these numbers refer to using the same ttg and the same 
parton distributions (fT3b irrespective of the order of the expansion. In data fits we expect a decrease 
of at NNLO of up to 10% with respect to the NLO value at about 30 GeV^ El El ESI . Also 
the parton distributions will be affected at this level, mostly in directions stabilizing the overall 
NNLO/NLO ratio Il24ll34ll35l . Thus such changes further support our conclusion of a rather good 
perturbative stability of Fi, for the jc- values accessible at HERA, at least above about 10 GeV^. 

Definite conclusions for F^ at considerably lower scales, say 2^ ~ 2 GeV^, require more de- 
tailed studies, since the larger value of ttg and the flatter small-x shape of the parton distributions 
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can lead to much larger NLO and NNLO corrections. This is roughly illustrated in Fig. 4, where 
we have employed (again independent of the order) the low-scale model distribution of ref. \2I, 



xg{x) 



0.6jc-«-1(1-x)3(1 + 10x'^-^) , 
1.2jc-0-ni--^)^l + 1.5x) , 



(14) 



together with tts = 0.35 for three flavours. Under these conditions the NNLO corrections now 
exceed 20% outside the range 7 ■ 10^^ < x < 0.1, rising sharply for decreasing lower values of x. 
Note that the low-scale gluon distribution is in fact poorly known at small x. Some of the recent 
analyses 041 l35l |9l |371 136| actually prefer a considerably smaller gluon density in this region, to 
the extend that at NLO can become almost zero, or even negative at 2^ < 2 GeV^ and very small 
X. More theoretical and experimental lfT2lfT3l efforts are required to clarify whether this behaviour 
and the large corrections in fig. 4 are analysis artefacts (in our case due to the unphysical order- 
independence of the inputs), or indeed indicate an anomalous low-order behaviour of Fi, or even 
signal an early (as compared to F2) breakdown of the perturbative expansion for this quantity. 

Fortran subroutines of our parametrized coefficient functions can be obtained from the 
preprint server http : / /arXiv . org by downloading the source of this article. They are also avail- 
able from the authors upon request. Routines for the exact results will be released with ref. ETl . 
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Figure 1: The perturbative expansion of the singlet-quark and gluon coefficient functions for Fi in 
electromagnetic DIS at a typical value of a^. The results have been divided by as = tts/ (47r), i.e., 
the LO curves are directly given by eq. ©. Also shown are the NNLO results at small x obtained 
if only the leading low-;c third-order terms |^|, xc\^ ^ ~ Inx, in eqs. Q and (fTOb are included. 
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Figure 2: The perturbative expansion of the singlet-quark and gluon contributions (up to a factor 
(e^), see eq. (HJl) to the longitudinal structure function at ~ 30 GeV^ for the typical (but 
order-independent) parton distributions (fT3t employed already for the illustrations in ref. Q- 
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Figure 4: As figure 2, but showing the results at a low scale ^2 GeV^ for accordingly modified 
values of ttj, njr and the (again order-independent) model distributions ([T4b used already in ref. [|2j|. 
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